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Abstract: We use recent results for the ^ Pl and — )• pT impact factors, computed 
in the impact parameter representation within the collinear factorization scheme, to get 
predictions for the polarized cross-sections ctt and 0"^ of the diffractive leptoproduction of 
the p meson at high energy. In this approach the helicity amplitude is a convolution of 
the scattering amplitude of a color dipole with a target, together with the virtual gamma 
wave function and with the first moments of the p meson wave function (in the transverse 
momentum space), given by the distribution amplitudes up to twist 3 for the 7^ — )• pT 
impact factor and up to twist 2 for the 72 — )• pL impact factor. Combining these results 
with recent dipole models fitted to DIS data, which include saturation effects, we show 
that the predictions are in good agreement with HERA data for photon virtuality (Q^) 
larger than typically 5 GeV^, without free parameters and with a weak dependence on the 
choice of the factorization scale, i.e. the shape of the DAs, for both longitudinally and 
transversely polarized p meson. For lower values of Q^, the inclusion of saturation effects 
is not enough to provide a good description of HERA data. We believe that it is a signal 
of a need for higher twist contributions in the p meson DAs. We also analyze the radial 
distributions of dipoles between the initial 7* and the final p meson states. 
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1. Introduction 

In the high energy hmit, exclusive processes and more particularly the diffr active lepto- 
production of vector mesons, provide a nice probe to study the hadronic properties. From 
the experimental side, data have been extracted in a wide range of center-of-mass energies, 
from a few GeV at JLab to hundreds of GeV at the HERA collider. The kinematical 
range which is at the heart of the present paper is the large energy in the center-of-mass 
of 7*p, denoted W, for which the HERA collaborations HI and ZEUS measured p-meson 
electroproduction starting from the early period of HERA activity [1,2] till now, with an 
increasing precision leading to a complete analysis [3, 4] of spin density matrix elements, 
polarized and total cross-sections describing the hard exclusive productions of the p and 
the (f) vector mesons V in the process 

j*{X^)p^ViXv)p. (1.1) 

These matrix elements and polarized cross-sections can be expressed in terms of helicity 
amplitudes T\y\^ (A^, Ay : polarizations of the virtual photon and of the vector meson). 

The ZEUS collaboration [3] has provided data for different photon virtualitics Q^, i.e. 
for 2 < Q2 < 160 GeV2, 32 < < 180 GeV < 1 GeV^), while the HI collaboration [4] 
has analyzed data in the range 2.5 < < 60 GeV^, 35 < < 180 GeV {\t\ < 3 GeV^) . 
The high virtuality of the exchanged photon allows the factorization of the amplitude into 
a hard subprocess described within the perturbative QCD approach and suitably defined 
hadronic objects, such as the dipole-nuclcon scattering amplitude or the vector meson 
wave functions and distribution amplitudes (DAs) [5-7]. HERA data are thus interesting 
observables to test the properties of these non-perturbative objects such as the saturation 
dynamics of the nucleon or the transverse momentum dependence of the vector meson wave 
functions. 

On the theoretical side, three main approaches have been developed. The first two, 
a A;T-factorization approach and a dipole approach, are applicable at high energy, W 3> 
Q S> A.QCD- They are both related to a Regge inspired fer-factorization scheme [8-14], 
which basically writes the scattering amplitude in terms of two impact factors: one, in 
our case, for the j* — p transition and the other one for the nucleon to nucleon tran- 
sition, with, at leading order, a two "Reggeized" gluon exchange in the i-channel. The 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution, known at leading order (LLx) [15-18] 
and next-to-lcading (NLLx) order [19-22], can then be applied to account for a specific 
large energy QCD resummation. The dipole approach is based on the formulation of sim- 
ilar ideas althought not in kx but in transverse coordinate space [23,24]; this scheme is 
especially suitable to account for nonlinear evolution and gluon saturation effects. The 
third approach, valid down to lU ~ Q, was initiated in [25] and [26]. It is based on the 
collinear QCD factorization scheme [27, 28]: the amplitude is given as a convolution of 
quark or gluon generalized parton distributions (GPDs) in the nucleon, the p-meson DA, 
and a perturbatively calculable hard scattering amplitude. GPD evolution equations resum 
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the collinear quark and gluon effects. The DAs are also subject to specific QCD evolution 
equations [29-31]. 

Though the collinear factorization approach allows us to calculate perturbative cor- 
rections to the leading twist longitudinal amplitude (see [32] for NLO), when dealing with 
transversely polarized vector mesons related to higher twist contributions, one faces end- 
point singularity problems. Consequently, this does not allow us to study polarization 
effects in diffractive p-meson electroproduction in a modcl-indcpcndcnt way within the 
collinear factorization approach. An improved collinear approximation scheme [33] has 
been proposed, which allows us to overcome end-point singularity problems, and which has 
been applied to /9-electroproduction [34-37]. 

In this study, we consider polarization effects for the reaction (1.1) in the high energy 
region, s = ^ ^ ^qcd^ working within the fcr-factorization approach, where 
the helicity amplitudes can be expressed^ as the convolution of the 7* — t- p impact factor 
^7*(At,)^p(Ap)^^2^ ^[^i^ ^i^Q unintegrated gluon density J^{x,k^) density, which at the 
Born order is simply related to the nucleon impact factor $'^~^^(fc^, Q^), where k is the 
transverse momentum of the t— channel exchanged gluons. For the transverse amplitude 
the end-point singularities are naturally regularized by the transverse momenta of the 
t— channel gluons [38-40]. At large photon virtuality, the 7* — p impact factors can be 
calculated in a model-independent way using QCD twist expansion in the region 
^QCD- Such a calculation involves the p- meson DAs as nonperturbative inputs. The 
calculation of the impact factors for (^Il-^pl^ (^It^pl jg standard at the twist-2 level [41], 
the next term of the expansion being of twist 4, while $Tt~^/'t was only recently computed 
[39, 40] (for the forward case t = tmin), up to twist-3, including two- and three-parton 
correlators, which contribute here on an equal footing. 

In a previous study [42], we used the results [40,41] for the ^"Il^pi- and (^'^t^pt 
impact factors and a phenomenological model [43] for the proton impact factor. It was 
pointed out that the region fc^ 3> ^^CD gives the dominant contribution to the helicity 
amplitudes while the soft gluon (fc^ < 1 GeV^) contribution cannot be neglected. The 
soft gluon contribution to the amplitudes involves the interaction of large size color dipole 
configuration jrj (|r| = l/|fc|) in the fluctuations of the probe and saturation effects could 
then play an important role. 

Following this idea, we have shown in rcf. [44] that the helicity amplitudes, expressed in 
impact parameter space and then computed in the collinear factorization scheme, factorize 
into the dipole cross-section and the wave functions of the virtual photon combined with 
the first moments of the p meson wave functions parameterized by the DAs, given by the 
twist expansion up to twist 3 for the production of a pT and to twist 2 when producing a 
PL. Note that in [38] a very similar approach was followed, where the DAs were replaced 
by the Taylor expansion of models of the p meson wave function at small qq pair transverse 
size. The results of ref. [44] link the ^t— factorization approach, in particular the results 
of [40], with the calculations performed in refs. [45-49] within the dipole approach. The 
main difference between our present approach and the one of refs. [45-49] is that instead 

^We use underlined letters for Euclidean two-dimensional transverse vectors. 



2 



of using the light-cone wave functions (f>{z,r) which in practice need to be modeled, the 
amplitude of our approach involves the DAs which parameterize the first moments of the 
wave functions. 

Our main point here is that one can assume the dominant physical mechanism for 
production of both longitudinal and transversely polarized mesons to be the scattering 
of small transverse-size quark-antiquark and quark-antiquark-gluon colorless states on the 
target. This allows to calculate corresponding helicity amplitudes in a model-independent 
way, using the natural light-cone QCD language - twist-2 and twist-3 DAs. 

This paper is organized as follows. In section 2, we introduce the impact factor repre- 
sentation of the helicity amplitudes as well as the kinematics of the process. In section 3, 
we first recall some results for the 7^ ^ Pl and 7^ px impact factors computed in mo- 
mentum space respectively up to twist 2 [41] and twist 3 [40] accuracies using the coUinear 
approximation to parameterize the soft part associated to the production of the p meson 
by distribution amplitudes (DAs), calculated in ref. [50,51]. We recall then the expression 
of these impact factors in the impact parameter space according to the results of ref. [44] 
allowing to decouple the dipole-nucleon scattering amplitude from the amplitude of pro- 
duction of dipoles in the initial (7*(A-y)) and final (p(Ap)) states. We terminate section 3 
by expressing helicity amplitudes and polarized cross-sections in terms of the dipole cross- 
section. In section 4, we give a brief review of the main properties of the models for the 
dipole cross-section [52-54] or the proton impact factor [43] that we use in our study. We 
compare our predictions with the data of HERA [3,4] in section 5 and we obtain a good 
agreement. In this context we discuss the role of higher twist corrections for small 
values. Finally, we analyse the radial distribution of dipole intermediate states involved 
between the virtual photon and the p meson, and discuss the role of the saturation models 
on the specific example of the Golec-Biernat and Wiisthoff saturation model [52]. We also 
compare our radial distribution with the overlap of the 7* and p meson wave functions 
obtained in the approach of the dipole models [46,47,55], where the p meson wave function 
is modeled and the parameters are fitted to HERA data. 

2. Helicity amplitudes of the hard p meson leptoproduction in the impact 
factor representation 



In the impact factor representation at the Born order, the amplitude of the exclusive process 
7*(A^) N p{\p) N reads 



as illustrated in figure 1. The 7*(A^) — )• p{\p) impact factor $t'('**7)-^p('^p) is defined through 
the discontinuity of the S matrix element for 7*(A-y; — >■ g{k — A)p{\p;pp) as 



Tx^x,{A;Q,M)=is 



I 



(fk 1 



-$^^^(fc,A;M2)$T*(^7)^^(^'')(fc,A;Q2)^ (2.1) 



$7*(A7)-^p(Ap) ^ ±_ 

2s 



P2 P2 ^ 



) 



(2.2) 
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where k = {k + q)^ . In eqs. (2.1) and (2.2) the momenta q and pp are parameterized via 
Sudakov decompositions, in terms of two hghthke vectors pi and p2 such that 2pi.p2 = s, 
as 

g = Pi p2 and p„ = H p2 + A_l , (2.3) 

s s 

where = —q^ » ^qcd virtuaUty of the photon being the large scale which 

justifies the use of perturbation theory, and nip is the mass of the p meson. Here — tmin 
denotes the minimal value of —t . The nucleon impact factor ^^^^ in cq. (2.1) cannot be 
computed within perturbation theory, and is related^ at Born order to the unintegrated 
gluon density J^{x,k). In the forward limit Aj^ = 0, the helicity amplitudes read, 
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/ '^-^^"y''{hQ,l^l)Hx.k). (2.4) 



The impact factors (^'i* ^^-^^^p^^p) Qz, Q, /x|.) and the nucleon impact factor vanish at ^ 
or — >■ A, which guarantees the convergence of the integral in eq. (2.4) on the lower limit^ 




Figure 1: Impact factor representation of the -y* N ^ pN scattering amplitude. 

The computation of the 7* — t- p impact factor is performed within collincar factoriza- 
tion of QCD. The dominant contribution corresponds to the 7^ — t- transition (twist 2), 
while the other transitions are power suppressed. The 7^ — )• Pl and 'j^ ^ pL impact 
factors were computed a long time ago [41], while a consistent treatment of the twist-3 
7^ pT impact factor has been performed only recently in ref. [40]. It is based on the 
light-cone collinear factorization (LCCF) beyond the leading twist, applied to the ampli- 
tudes j*(Xy)g(k) —7- g{k — A)p(Ap), symbolically illustrated in figure 2. Each of these 
scattering amplitudes is the sum of the convolutions of a hard part (denoted by H and Hp 
for two- and three-parton contributions, respectively) that corresponds to the transition of 
the virtual photon into the constituents of the p meson and their interactions with off-shell 
gluons of the t channel, and a soft part (denoted by $ and $'*) describing the hadronization 
of the constituent partons into the p meson and parameterized up to a given twist by DAs. 



^Normalization of the impact factors differs from [56] by a factor 2n, $ = 27r^^'^^''. For clarity, we 

e restored in eq. (2.4) the colore 

^This property of the impact fac 
related to gauge invariance [57,58]. 



have restored in eq. (2.4) the colored indiced carried by the impact factor 3>^j . 

^This property of the impact factors is universal in the case of the scattering of colorless objects and is 
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Figure 2: Two- and three-parton correlators attached to a hard scattering amphtude in the specific 
case of the 7* — >■ p impact factor, where vertical lines are hard t— channel gluons in the color singlet 
state. 

3. Helicity amplitudes and polarized cross-sections 
3.1 Impact factors 7^ 7- — > pl,t 

In the Sudakov basis, the longitudinal and transverse polarizations of the photon are^ 

<L = ^(Pi+ V^2) , = ^ (0, Tl, -i, 0) . (3.1) 

For t = tinin the same parametrization will be used for the /9-meson polarization with 
obvious substitutions — t- — and Q — t- nip. 

Let us introduce p and n, two light-cone vectors such that pp ^ p at twist 3 and 
p - n = 1. The polarization of the out-going p meson is denoted by e*. For further use, we 
denote R]^^ = Saxps e*l^p^ ■ We use the following convention for the transverse euclidean 
polarization vectors as in eq. (3.1) 

e± = -^(Tl,-z) . (3.2) 

The DAs needed to parameterize the p meson productions involved in the results of [40] 
for the impact factors with A = 0, according to [40] are: the twist 2 DA (pi{y;p'jp) associ- 
ated to the production of the longitudinal p meson, the two-parton twist 3 DAs ^Jiy, fJ-'p), 
(p^{y; p?p) and the three-parton twist 3 DAs B{yi,y2; p'jp) and D{yi,y2; Pp) that parame- 
terize the production of a transversely polarized p meson. Following [42,44], we also use 
the combinations, 

M{yi,y2;PF) = (pp) B{yi,y2; pp) - CiiPp) D{yi,y2; Pp) , 
S{yi,y2\pl) = CUpl)B{yi,y2;pl) + C3{pl)D{yi,y2;pl), (3.3) 

where CXit^^^) Czit^V) the dimensionless coupling constants 

(3^1) = :^, (3^(4) = ^%^. (3.4) 

The scale pp is the factorization scale involved in the production of the p meson, that we put 
equal to the renormalization scale of the evolution of the DAs. We recall in appendix A.l 

''in ref. [40] we took = (0, 1, ±i, 0), which we change here for consistency with the usual experi- 
mental conventions [59]. 
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some basics features of the ehiral even twist 2 and twist 3 DAs present in this approach; 
more details can be found on the DAs in ref. [40]. We recall also in appendix A. 2 the 
dependence of the DAs on the collinear factorization scale /ip which is driven by the 
renormalization evolution equations given in ref. [50] . 

We have recently shown in [44] that these impact factors, expressed in the impact 
parameter space, read, 

^il^PL = -Impfpj dy J d\{e* ■n)ipi{y-iJil)Hll{y,r,k), (3.5) 

and 

^7|-^PT = -^mpfpj dy2 J d\{ip3{y2;p,l)Ht}{y2,r,k) + iipJ{y2-,,J.l)ie*-r)H^^^^ 
+ i'PAiy2;f^l)H^^^'iy2,r,k) - vl{y2\pil) (E* ■ l) H^gf'{y2,r,k) 

+zC3^i?(yi,y2;/x|^)^S''*''''(yi,y2,r,r',fe))} , (3.6) 

where we have respectively denoted in the case of two-parton exchange as y and y = \ — y 
the longitudinal momentum fractions of the quark and the antiquark, and in the three- 
parton exchange yi , ^2 = 1 — 2/2 and yg = y2 — yi the longitudinal momentum fractions of 
the quark, the antiquark and the gluon. The transverse displacement vectors of the color 
dipolc configurations are denoted as r for the interacting dipole (two- and three-parton 
Fock component) and r' for the spectator dipole (present only in the three-parton Fock 
component). Note that jrj in the case of the two-parton component is the transverse size 
of the quark anti-quark colorless pair. The functions denoted 

m ~ m 1^ ' 999 ~ 999 " I" ' J 

are the Fourier transforms in the transverse plane of the two-parton component hard parts 
H and the three-parton component hard parts H^j, (illustrated in figure 2), projected on 
the appropriate Fierz structures P^. 

The computations of the hard parts lead to the following generic expressions, 

^^l^P^ {k,Q,^il)= [^-^^ j dy j dr V;^^^''^ (y, r ; Q, /.|) A{r, k) , (3.8) 

^^T^P^k,Q,txl) = j dy j dr ^^J^'"(y,r;Q,/x|,)^(r,^) (3.9) 

( V) / / '^^^ j '^^'^S^''^(yi'^2,2:;Q,/^F).A(E,:fc), 

where the function 

A.'jrf)/ 

Air.k) = (1 - exp {ik ■ r)) (1 - exp {-ik ■ r)) , (3.10) 
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is the scattering amplitude of a quark- antiquark color dipole with the two t— channel glu- 
ons, putting apart the color factor Tr(t" t^) = (^"^/2, with a and b color indices and Nc the 
number of colors. The functions '>pjgq^^^ j '^Jqg)^^'^ ' '^'{qqg)'^ ^'■^ respectively the amplitudes 
of production of a p meson from a quark-antiquark (quark-antiquark gluon) system pro- 
duced far upstream the target in the fluctuation of the virtual photon. These functions are 
computed up to twist 3 in the collinear aproximation in ref. [44] and they contain informa- 
tion about the relevant color dipole system that interacts with the target. The functions 
^(s^^r expressed in terms of the virtual photon wave functions ^^j^'f^y 

^jyy,r;Q') = h,-H^^^K,{,.\r\) , (3.11) 

*Ji)^^(y,r;Q^) = h,-H'i^^{yh,x, + y6h,-x,) ' lul) , (3.12) 

where h = h = zt^ denote respectively the helicities of the exchanged quark and 
anti-quark, and of the combinations of DAs of the p meson ^1*^^^ , 



= h-h\j^{(^L ■ , (3-13) 



kPL I,.. ,,2 
\h,h) 



X i4) + {Sh,x, - Sh,-xMiy' ^'F)) ■ (3-14) 

In the two-parton approximation, these functions (^^^^^(y; /Xp) and 4>1j^'^^''\y,r; /j,'^) pa- 
rameterize the moments of the wave functions of the p meson, i.e. the first terms of the 
Taylor expansion of the wave functions at small r. Note that in the approach of rcf. [38] 
'^{hh)^y'^ ^"f^ ^^'^ '^{hh)'''^^y^—- ^'''^f) replaced by the Taylor expansion for small r of the 
modeled wave functions of the vector mesons. Finally, the functions ''pjqq^'^^ ^(99)^^^ 
with two-parton components read 

4^'Hy,r;Q,^^l) = ^ E <.%(^;Ml)*5,,)(y,z:;Q^), (3.15) 

(h,h) 

^t,7^{y,r,Q,,l) = ^ E ^(^;Ml)*Ji;^^(y,i:;Q^). (3.16) 

{h,h) 

The function i^J^q^^^^ with three-parton components reads 

,2 



^(;,,) (y,r;Q,/x^) = ^(^^— ^ J ^-(?/i,y2,r;Q) 



J^^{y2,yi,r; Q) 



(3.17) 



where the function T'^t describes the fluctuation of the transversely polarized photon into 
a quark-antiquark-gluon color singlet. The function J^'^t can be expressed in terms of the 
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longitudinally polarized photon wave function 



{h,h) 



27rV TT Q 



-Kq^ih \r\ 



(3.18) 



as 



^^(yi, 2/2,1; Q) 



viQ 



+ 



+ 



Nc 
Cf 
Nc 
Cf 



yiQ 



+ 



- 2 



+ 



y2yij \ yiQ 

ygQ ygQ 
y2 f^''Hf^2,r;Q) 



yiQ 



y2 



ygQ 



ygQ 



with 



Ml = yiyiQ^ , M2 = y2y2Q^ 

,9 _ ViVg ^2 „2 _ _mi^^2 

yi + yg y2 + yg 



(3.19) 



(3.20) 



yiy2 ^2 



yi + y2 

and Cf = 2N ■ Note, that in the large Nc limit 7-"^? simplifies, 

1 



J^^iyi,y2,r; Q) 



X {yiy2*^MMi,z:;Q) + y2j/i*^^(/X2,z:;g) -yg^'^i(/Xgg,r;Q)} . 

3.2 From impact factors to helicity amplitudes and polarized cross-sections 

According to eq. (2.4), the helicity amplitudes read, 



(3.21) 



(3.22) 



Inserting the expressions for the impact factor i^'^^i^'^^p of eqs. (3.8, 3.10), one gets 

~J' ^ J J dr^pjqg)'''^ (y^l-'^Q^ f^p) ^{x,r) 
where we have identified the dipole cross-section as defined in ref. [56] 



(3.23) 
(3.24) 



(T(x,r) = 



N^-1 f cfk 



(3.25) 
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Note that we can separate the Tn as the WW contribution and the genuine contribution, 
^ = jdyj dr^J^^''-''^'^{y,r,Q,^^l)aix,r), (3.26) 

rpgen f f * 

-^ = JdyJ dri;J-^'-'''''{y,r,Q,^^l)aix,r) (3.27) 
+ Jdy2j dyi J dr%l)J^^^'''' {yuy2,r,Q, i?p) a{x,r) , 

The formulas (3.23, 3.24) allow us to combine various models of the scattering amplitude 
of a dipole on a nucleon with the results obtained by twist expansion of the 7* — > p impact 
factor. At t = tmin the contribution to the cross-sections ctl and gt are respectively coming 
from the helicity amplitudes Too and Tn, 

^,^,^EW^. ,,,, 

The t— dependency is expected to be governed by non-perturbative effects of the nucleon, 
which can be phenomenologically parameterized by an exponential dependence of the dif- 
ferential cross-section 

^(i) = e-^W^)*^(i = 0). (3.30) 

This results in the polarized cross-sections 

1 |Too(^,t = 0)|^ 
^"^ = 6(0^ 16^^ ' ^^-^^^ 

1 |rn(^,t = o)|^ 
^^ = 6(g^ — 16^^ — • ^^-^^^ 

The b{Q'^) slope has been measured by ZEUS and HI. We will use here quadratic fits of 
the b{Q^) slope data of ref. [4] to determine the cross-section. 

In the following section, we will briefly present the dipole models we shall use to 
compare our predictions with HERA data. Note that we could in principle use a dipole 
model taking into account skewness effects. This skewness dependence can be implemented 
along approaches of refs. [60,61], but this subleading physical effects will be neglected in 
the present study. 



4. Dipole models 

In the dipole picture, the DIS cross-section reads [23, 24] 

r r I * 2 

<T = J dyY,\^'riy,L;Q) Hx,r), (4.1) 
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with 



a{x,r_) = j dH^ = 2j d%N{x,r,h). 



(4.2) 



The low-x saturation dynamics of the nucleon target was first introduced in refs. [52, 62] 
by Golec-Biernat and Wiisthoff (GBW) model to describe the inclusive and diffractive 
structure functions of DIS, which inspired many phenomenological descriptions of DIS 
HERA data [47,55,63-67]. In this model the dipole cross-section reads 



a(x, r) = o"o < 1 — exp 



4i?^(x) 



where 



Rl{x) 



1 



GeV^ V.XO 



X 



(4.3) 



(4.4) 



and it involves three independent parameters {(To,xo, Ap}. One can see from eq. (4.1) that 
since the wave functions are peaked at r ~ ^, the domain in which the saturation effects 
are significant is given by 

1 o2, 



(4.5) 



To make contact with photoproduction, it is customary [52] to make the following modifi- 
cation of the definition of the Bjorken variable x 



X — )• X 1 + 



1 + 



Q2 I Q2^o ' 



(4.6) 



where m/ is an effective quark mass which depends on the flavour / and of the model used 
to fit the data. This modification is adopted in the following. 

In the GBW saturation model, light quark masses are taken to be 0.14 GeV in order 
to get a good fit of the photoproduction region. The inclusion of the charm contribution, 
with rric = 1.5 GeV has also been performed in ref. [52]. The normalisation ctq results from 
the integration over the impact parameter b, assuming that the 6— dependence of the dipole 
amplitude factorises as M{x,r,b) = T{b)J\f{x,r), where T{b) is related to the density of 
gluon within the target, leading to 



a{x,r) = 2 (fbT{b)M{x,r) = aoM{x,r) 



(4.7) 



This model provides a good description of inclusive and diffractive structure functions for 
the values of the parameters presented in table 1. 



Fits 


cTo (mb) 


Xp 


Xq 


No charm 


23.03 


0.288 


3.04 X 10-^ 


With charm 


29.12 


0.277 


0.41 X 10-^ 



Table 1: Values of the parameters entering the GBW dipole cross-section. 
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The small x evolution of the dipole cross-section can be modeled as in the refs. [52], 
or computed numerically from the running coupling Balitsky-Kovchegov (rcBK) equa- 
tion. Indeed, the a;— dependence is driven at small x by perturbative non-linear equa- 
tions, the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) equation 
[68-73] and the Balitsky-Kovchegov (BK) equation [74, 75] . The solutions of the two equa- 
tions are not significantly different and the BK equation is simpler to solve numerically 
than the JIMWLK equation. The LO-BK equation cannot be used to obtain the low x evo- 
lution as it predicts a growth of the saturation scale much faster than the growth expected 
from the anlysis based on phenomenological models. It was shown in [53, 54] that taking 
into account only the running coupling corrections of the evolution kernel allows to get the 
main higher order contributions, and to solve the discrepancy between the growths of the 
saturation scale. A numerical solution of the running coupling BK (rcBK) equation was 
obtained in refs. [76, 77] based on initial conditions inspired by the GBW model [52] and 
the McLerran-Venugopalan (MV) model [78]. We will denote these numerical solutions 
for the dipole scattering amplitudes as the Albacete-Armesto-Milhano-Quiroga-Salgado 
(AAMQS)-model. 



AA*^V'(xo,r) = a^^^|l-exp 



4 

2/n2 \ 7 



In 



Aqcd 



+ e 



(4.8) 
(4.9) 



with xq = 0.01 and Qso the initial saturation scale at x = xq, and 7 the anomalous 
dimension. The coupling constant in the evolution kernel of the rcBK equation depends 
on the number of active quark flavors Uf, 



-Itt 



4C2 



(4.10) 



where Po,nf = H ~ l*^/) is the QCD scale and C is one of the free parameters of the 



model. As usual, the scales A„^ are determined by the matching condition 



Q^s.rif — 1 



(ri) 



"s,n/(?'*) at rl = iC /mj: and an experimental value of a^. In eq. (4.9), the scale Aqcd is 
identified with A3. 

The parameters are fitted to the experimental data for the inclusive structure function 
of DIS 



F2{X,Q') 



47r2 



Op 



and give a good description for the data of the longitudinal structure function 



Fl{x,Q' 



47r2 



-CTL ■ 



(4.11) 



(4.12) 



Two sets of dipole cross-sections are available for each of the initial conditions. The first 
set of dipole cross-sections respectively denoted as (a) and (e) for the GBW and the MV 
initial conditions, are fitted to data by taking only into account the light quarks u,d,s. The 
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second set, respectively denoted as (b) and (f) for the GBW and the MV initial conditions, 
are fitted taking into account the dipole cross-section of the light quarks aoAf^'{x,r) and 
of the heavy quarks aQj\f^'{x,r) . The normalization of the dipole cross-section erg' for the 
charm and the bottom quarks are assumed to be equal. Thus the cross-sections read 



a 



.TP 



L,T; set(a),(e) 



,7 P 



L,T,set(b),(f) 



f=u,d,s 

= 4 [d\[dy\^f'^iy,r;Q,mf,ef)fAf\x,r) 

f=u,d,s 

+ ^0 I^^'^J dy\^f'^{y,r;Q,mf,ef)f M''-{x,r). 



(4.13) 



f=c,b 



(4.14) 



In the following parts we will focus mostly on the sets (a) and (b) using the GBW initial 
condition respectively denoted as AAMQSa and AAMQSb, as the effects due to different 
initial conditions do not involve significant changes in the numerical results of present 
study. We present in Tabs. 2 and 3 values of the parameters of the fits obtained in ref. [77]. 



Fits 




(To (mb) 


7 


C 


xVNaf 


(a) 


0.241 


32.357 


0.971 


2.46 


1.226 


(e) 


0.165 


32.895 


1.135 


2.52 


1.171 



Table 2: Values of the parameters entering the AAMQS sets (a) and (e) dipole cross-sections. 



Fits 


Q% 


M) 2 


o-Q (mb) 


a^- (mb) 


7 




C 




(b) 


0.2386 


0.2329 


35.465 


18.430 


1.263 


0.883 


3.902 


1.231 


(f) 


0.1687 


0.1417 


35.449 


19.066 


1.369 


1.035 


4.079 


1.244 



Table 3: Values of the parameters entering the AAMQS sets (b) and (f) dipole cross-sections. 



Note that the models of dipole cross-section we use here involve only a quark anti-quark 
intermediate state between the initial and final wave functions used to describe DIS process 
while the impact factor — > pr involves also qqg partonic state. The error induced by 
this approximation should be subleading compare to higher twist corrections or the choice 
of fiF, but still it is difficult to appreciate quantitatively how much it could impact on the 
predictions. 

For completeness, we will also display predictions using the Gunion and Soper model 
(GS-model) [43]. This model was used in our first phenomenological study of ref. [42], and 
it assumes that the hadron impact factor takes the form 



1 



M2 + (|)2 M2 + (^-|)2 



(4.15) 



where A and M are free parameters that correspond to the soft scales of the proton-proton 
impact factor. As it was discussed in [42], the ratios of helicity amplitudes are well described 



12 



for M ~ 1 GeV and the result is not very sensitive to this parameter. Note that this impact 
factor indeed vanishes when fe->OorA — fc— )'Oina minimal way. Such a model was 
the basis of the dipole approach of high energy scattering [79] and used successfully for 
describing deep inelastic scattering (DIS) at small x [80]. The dipole-proton scattering 
amplitude in the two-gluon exchange approximation, computed in details in appendix A.3, 
resulting from the GS-model for the proton impact factor is: 

= / # (FjJ -^*«' {w - M^) (' - ''"^) (' - 

5. Compctrison with the HERA data 

Our main results are the polarized cross-sections ut and uj^ of the process (1.1), that we 
compare with the data of the HI collaboration [4]. 

In the following plots, experimental errors are taken to be the quadratic sum of sta- 
tistical and systematical errors. The coUinear factorization scale jip only appears in the 
scattering amplitudes through the DAs and the coupling constants and unless specified we 
will assume that jip depends on the virtuality as 

,|(Q=,.«!±^^?^^^. (5.1) 

The figures 3(a) and 3(b) show the full twist 3 predictions respectively for gt and 
(Tl, using the different dipole models. The results of the predictions for ctt (figure 3(a)) 
and (Tl (figure 3(b)) are in agreement with the data for values of larger than approxi- 
mately 7 GeV^ for cry and 7 GeV^ for at depending on the considered dipole model. The 
results based on the AAMQSa model are giving a slightly better description than with 
the AAMQSb although both sets give good quality fits of DIS data. The results obtained 
by GBW model are close to the one obtained by the AAMQSa model; both are in good 
agreement with the data for above 7 GeV^. Let us insist on the fact that the agreement 
for > 7 GeV^ of our predictions with data is non-trivial since all free parameters of the 
dipole cross-section (normalization, Rq{x), etc..) are completely fixed by the fit of DIS 
data, while on the p side the normalizations are given by decay constants obtained from 
QCD sum rules. 

Keeping in mind that the results are not too sensitive to the precise choice of the dipole 
model, below we will focus on the predictions of the AAMQSa model. Some results of the 
GBW, AAMQSb and GS models are presented in appendix A. 4. 

In figure 4, we show separately three different contributions: 

• the full twist 3 (Total) contribution, involving both the WW and the genuine solutions 
of the DAs. 

• the WW contribution, only involving the WW solutions of the DAs. 
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cr J (nb) 




(a) (tt vs HI data 



cr^(nb) 




(b) (TL VS HI data 



Figure 3: Predictions for ctt and (7^ vs Q^, for W = 75 GeV, using the AAMQSa (red solid line), 
AAMQSb (blue large dashed line) and GBW (green dashed line) models compared to the data of 
HI [4]. 

• the asymptotic (AS) contribution (for fip — t- oo), involving the asymptotic solutions 
of the DAs. In this limit the genuine contribution vanishes and the WW DAs can be 
expressed as functions of the asymptotic DA (pi (y) = 6yy. 

The results of the predictions for ctt (figure 4(a)) and (Jl (figure 4(b)) are in agree- 
ment with the data for values of respectively larger than Q'^"^^^ ~ 6.5 GeV^ and 
Q2min ^ g GeV^ which confirms that the amplitude factorizes into a universal colour 
dipole scattering amplitude and that the truncated twist expansion of the p meson soft 
part is justified. 

The two scales Q^™'° and Q^™™ are close to each other. We interpret this fact as an 
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cr J (nb) 




(a) AS (purple dashed line), WW (blue long dashed line) and Total (red solid line) 
contributions to (Tt- 



o-L(nb) 




(b) AS (purple dashed line) and Total (red solid line) contributions to ctl. 



Figure 4: Predictions for (Tt and cr^ vs Q^, for W — 75 GeV, using the AAMQSa-model, compared 
to the data of HI [4]. 

indication that the discrepancy between data and our predictions at low are mainly 
due to the higher twist contributions to the impact factor from the meson structure rather 
than an effect of the saturation dynamics of the nucleon which should be well described at 
this scale by the saturation models. Let us mention, again that these saturation models 
are known to fit very well inclusive DIS as well as diffractive DIS data (for GBW) at these 
low values. 

The saturation scale, given by Qs = l/Ro{x) is of order 1 GeV in the kinematics of 
HERA. Since our predictions are only consistent with data in the region > Q|^™™ > Q|, 
this limitation do not allow us to access the domain Q| > where saturation effects can 
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be essential. 

The predictions are dominated by the WW-contribution and are not very sensitive to 
the choice of the coUinear factorization scale. This will be further discussed in section 6. 




Figure 5: Quadratic fits of the 6— slope HI data. 



An estimation of the error on the cross-sections caused by the error bars on the 6— slope 
measurements is obtained by fitting the upper and lower bounds of the 6— slope as shown 
in figure 5, and then by computing the predictions based on these fits as shown in figure 6. 
Note that we have assumed that the longitudinal bL and the transverse bx slopes are equal 
to the 6— slope of the total cross-section. This assumption is supported by HI data where 
the measurements of the difference bi — bx for = 3.3 GeV^ and = 8.6 GeV^ are much 
smaller than the 6— slope value. Let us also emphasize that in this approach we compute 
the polarized differential cross-sections in the limit t = tmin ~ where only the s-channel 
helicity conserving (SCHC) amplitudes Too and Tn are non-zero. The contributions of other 
helicity amplitudes are encoded in the phenomenological t— dependence given in eq. (3.30) 
and it turns out that data for the total differential cross-section are dominated by a t-region 
of very small values, with a typical spread given by the scale ^ \ GeV^. 

We now compare our predictions with the data for the total cross-section a, given by 
the sum a = ai + c^t according to ZEUS convention in ref. [3] ov a = eai + cft following 
HI notation [4], where e is the photon polarization parameter^ 

e~ (l-y)/(l-2/ + yV2). (5.2) 

We show in figures 7(a) and 7(b), the AS, WW and Total contributions to the total cross- 
section a as function of for fixed averaged W . The predictions are larger than the 
data for smaller than approximately 7 GeV^, as expected from the results of gl^t- The 
figures 8(a), 8(b), 9(a) and 9(b) show the W dependence of the cross-section for several 
values of Q^. We see again a good agreement between the predictions and the data for 

^Yor HI (e) = 0.98 and for ZEUS (e) = 0.996. 
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cr J (nb) 




(a) AS (purple) and Total (red) contributions to gt- 



cr^(nb) 




(b) AS (purple) and Total (red) contributions to ctl. 



Figure 6: Full twist 3 and asymptotic predictions with the 6— slope uncertainty, using AAMQSa 
model. 



approximately above 6 GeV^ for HI data and 8 GeV^ for ZEUS data, taking into account 
the uncertainty on the 6— slope. 

Finally, our analysis also provides predictions for the ratios R and for the spin density 
matrix element Tqq 

R=—, (5.3) 
(It 

r'o'o = ^ • (5.4) 
Assuming to keep only the SCHC amplitudes Tn and Tqq, and the equality of slopes 



17 



o-(nb) 




(a) AS (purple dashed line), WW (blue long dashed line) and Total (red solid line) 
contributions vs HI data. 



(T(nb) 




(b) AS (purple dashed line), WW (blue long dashed line) and Total (red solid line) 
contributions vs ZEUS data. 



Figure 7: Predictions for a vs compared respectively with HI [4] data (figure (a)) for W — 
75 GeV and with ZEUS [3] data (figure (b)) for = 90 GeV, using the AAMQSa-model. 



= ^T, the i— dependences of the cross-sections cancel in the ratios R and Tqq, leading to 



R = — (5.5) 

^11 



and 



^00 = , (5.6) 



where xn = |Tii|/|Too| . HI and ZEUS measurements of R and Tqq = ctl/ct as functions 
of |t| confirm this weak dependence on \t\. Based on HI data, we can estimate [42] the 
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Q- = 6.6 GeV^ 



o-(nb) 

1000 f 
500 



100 
50 

10 

5F 



2" = 35.6GeV^ 
W(GeV) 

60 80 100 120 140 160 180 

(a) Asymptotic (purple dashed line), WW (blue long dashed line), total (red solid 
line) contributions vs HI data. 



= 3.3 GeV^ 




40 



= 6.6 GeV^ 




W(GeV) 

40 60 80 100 120 140 160 180 

(b) Total contribution to a including the fe— slope errors vs HI data. 

Figure 8: Predictions for cr vs W, using the AAMQSa-model, compared with HI [4] data. Figure 
(a): AS, WW and Total contributions. Figure (b): Total contribution taking into account the 
uncertainties on the 6— slope. 

correction to the ratio Tqq due to the amphtude Tqi, for the range of HI, to be below 
1%. The results are shown in figure 10 for the ratio R and in figure 11 for the spin density 
matrix element Tqq, using AAMQSa model. 



6. The radial distributions of dipoles involved in the overlap of the 72(r) 
and pi(r)— meson states 



In the dipole picture, the overlap of the wave functions of the outgoing p meson ^' 



*Pl,t (Ap) 
(h,h) 
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(r(nb) 
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^ = 3.7 GeV^ 
s* Q- = 6. GeV- 




= 13.5 GeV' 



j2 + 



1 1 2 - 32. GeV" 

^ ^ ^ ^ — ' ' ' ' ' ' ' ' ■ ' ■ ' W(GeV) 

30 50 70 100 150 200 

(a) Asymptotic (purple dashed line), WW (blue long dashed line), total (red solid 
line) contributions vs ZEUS data. 
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(b) Total contribution to a including the 6— slope errors vs ZEUS data. 

Figure 9: Predictions for a vs W, using the AAMQSa-model, compared with ZEUS [3] data. 
Figure (a): AS, WW and Total contributions. Figure (b): Total contribution taking into account 
the uncertainties on the 6— slope. 



and of the incoming virtual photon ^'Jf^'?/ represents the amplitude of probability for 
these states to dissociate into a quark anti-quark color dipole of size r, the quark having 
a longitudinal momentum fraction y. We define thus the probability amplitude Wa^a^ as 
the corresponding parts of the impact factors appearing in eqs. (3.23) and (3.24), 



(6.1) 



20 



O-JO-T 

10 r 



Q- = 22.5 GeV- 

Q- = 7.5 GeV"_ i 
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(a) AAMQSa versus HI data. 



W(GeV) 



Q- = 18.0GeV 

g- = 8.8 GeV- 
2^ = 4.2GeV 
= 2.4GeVi 




50 100 150 

(b) AAMQSa versus ZEUS data. 



W(GeV) 



Figure 10: The full twist 3 contribution to the ratio of the cross-sections R — UhjcrT in the limit 
i = versus W and compared to the data of HI [4] in figure (a) and ZEUS [3] in figure (b). 

for the amplitude Tqo and 



^(gg) \y ■,!-■, ^ ■• F ) ■• (6-3) 

,^7j, ^PT gen ^.ri ,,2^ 1 / " ^a.. jCt^Pt ( „, „, ^. ,.2 

/o 

for the Total, the WW and the genuine contributions to Tu. The probability amplitudes 
WapAt permit in turn to define the radial distributions Vx^x^ of the interacting dipole, as 



7^A,A,(r-,Q',MF) = T7 r / dy \Wx,x,{y,r- ^il,Q^ , 



where N'x„x^, are normalization factors 



AAa 



drr / dy \Wx„xJyi,y,r; fiF,Q ) 



Expressed in terms of these functions, the scattering amplitudes read 



AAa 



POO 

/ drVxpX^{r,Q'^,fj.p)a{x,r) 
Jo 



(6.5) 



(6.6) 



(6.7) 



The probability amplitude for a dipole of size r to scatter on the nucleon is then proportional 
to VxpX-yi^^Q'^ T IJ-f) ^{x,r), which justifies the inclusion of the factor r in (6.5). 
Below we will also use the rescaled radial distributions Px^x^i^^ fJ'F)^ 



PXpXJX,f^F)^ 



2, 'PxpX^{^,Q'^; fJ-p) 



Q 



(6.8) 



which depend on r and Q only through the variable X = rQ and we choose to put //^ 
//p((5^), see eq. (5.1). Note that the rescaled asymptotic distributions. 



(6.9) 
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(a) Asymptotic (purple dashed line), WW (blue long dashed 
line), Total (red solid line) contributions using AAMQSa-model 
vs HI data. 
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(b) Asymptotic (purple dashed line), WW (blue long dashed 
line). Total (red solid line) contributions using AAMQSa-model 
vs ZEUS data. 

Figure 11: Predictions for Tqq vs W and compared respectively with HI [4] data (figures (a)) 
and ZEUS [3] data (figures (b)), using tire AAMQSa-model. 



are independent of Q^. 

This change of variable leads to the formulas 



poo 

^fx,x, / dAPx,A,(A,//|)a(x,A), (6.10) 
Jo 

with 

a{x,\)=a(^x,^^ . (6.11) 

The average value of a function /(y) depending of the longitudinal fraction of momen- 
tum y carried by one of the partons, will be estimated by 

if(y))x,x, = fdrf dyf{y)r \Wx,x,{y,r; i^l,Q^)\ . (6.12) 
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0.1 0.2 0.5 1.0 2.0 5.0 10.0 0.1 0.2 0.5 1.0 2.0 5.0 10.0 

(a) Total Poo(r,Q^A^F(Q^)) (red solid), (b) Total Poo(r, Q^ /iF(Q^)) (red solid), 

Vl)o'''\r,Q'^) (red dashed) distribut ions and 'Pqo^\i~,Q^) (red dashed) distributions and 

a{x,r) at = 1 GeV^ a{x,r) at = 10 GeV^ 



Figure 12: The functions Voo[t,Q^ , ^Pp) and T'qq (r, Q^) vs the size r of the interacting dipole, 
for Q'^ — \ GeV^ (12(a)) and = 10 GeV^ (12(b)), and the dipole cross-section (7(x,r) normahzed 
by the factor 4cro for = 50 GeV and W = 150 GeV. 



6.1 The radial distribution of the 72 pi transition 

The distributions Vooir, Q'^^fil) and V^oQ^\r, Q^) = Voo{r, , oo) are close to each other, 
as it is shown in figures 12(a) and 12(b), which indicates that the distribution Voo{r, Q'^, fip) 
is not sensitive to fip. We can then restrict the study of Voo{r, Q"^, fip) by considering only 
'PQQ^\r,Q'^), which is also simpler for the analytic treatment. At first glance we see that 
the distributions are peaked around ~ ^ and consequently the peak moves to the right 
and the distribution becomes wider as decreases. Note that the dependency of a{x,r) 
with respect to Q^, which can be seen in figure 12, only occurs through the dependency of 
Ro{x), according to eq. (4.4). 




0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 



Figure 13: Full twist 3 (Total) Pqq{X, ^ll,{Q^)) for = IGeV^ (solid red) and = lOGeV^ 
(dashed blue), AS P^o'^\x) (dotted purple) and a{x, A) at 1^ = 90 GeV^ for = 1 GeV2(dotted- 
dashed black) and = 10 GeV^ (dashed black). 

In figure 13 we show the Total and AS rescaled radial distributions Pooi^, l^'piQ'^)) ^-nd 
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P(S^'^V>Mf(<3^))- This last one reads 

Pif\>^) = ^'Pl,f\^,Q') = 6 I dy{yyf\Ko{VW>^). (6.13) 

The average value of A estimated with Pqq^\\) is 

(A)if)=/dAAPo^^)(A) = ^.3.7. (6.14) 

About half of the dipoles are contained in the region 1 < A < (A)q^'^'' , the peak of the 
distribution being at AP^^'' ~ 1.3. The typical transverse scale /x = \/yy entering the 
wave functions overlap can be estimated using eq. (6.12), 



{AS) 
00 



Q 



= <V^>or^ =^ l^^y I ^^(V^) {{yyfXKo{VyV>^)) = 0.44. (6.15) 



The choice of the factorization scale ^ipiQ ) given by eq. (5.1) is then a good approximation 
of the transverse dynamical scale (/^)oo'^^ involved in the process. 

The dipole scattering amplitude plays the role of a filter that selects dipoles of A > 
A^^* (Q^,PF) = 2Ro{x)Q. Note that the critical saturation line eq. (4.5) is given by 

\''^'-{Ql{x{Q%W^)),W) = 2, (6.16) 

where 

in accordance with eq. 4.6. In the kinematics of HERA, the energy in the center of mass 
W varies roughly from 50 GeV to 150 GeV, leading to the following bounds for the two 
values g2 = 1 GeV^ and = 10 GeV^, 

A^^*-(l,50) = 2.1 > X^^^-{1,W) > A^'^'- (1,150) = 1.5 (6.18) 
A^^* (10,50) = 9.2 > X^^^-{W,W) > A^^*- (10, 150) =6.7. (6.19) 

We will fix for our purpose W = 90 GeV resulting in the values, (AS^*-(1 GeV^90 GcV)~ 
1.8) and (AS^* (10 GeV^ 90 GeV)~ 7.7). We can then differentiate the case = i GeV^ 
where we are in the saturation regime (A^^* (l GeV2,90 GeV)< 2), and the case = 
10 GeV^ where saturation effect are less important (A^^* (10 GeV2,90 GeV)> 2). 

We can evaluate the percentages Nx^xsi^t. of the dipoles large enough to be in the 
bandwidth of the dipole cross-section, for each Q"^ , 



/•OD 

Nx>xs-t- {Q^ = 1 GeV^ = 90 GeV) = / dX P^^^' (A) = 70% , (6.20) 

7ASat(l,90) 

POO 

Nx>xs.t. {Q^ = 10 GeV^ W = 90 GeV) = / dX P^^^^ (A) = 10% , (6.21) 

JASat (l0, 90) 

as one can see in figure 13 (plotted in logarithmic scale). The large difference between 
-^A>ASat. (1, 90) and A^;)^>;)^sat. (10, 90) indicates that the integrand of Too shown in figure 14, 
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Figure 14: The normalized integrand of Too, i-e. Voo{r,Q^ , fi^p) a{x,r). The Total integrand at 
/i|,(g2) for Q2 = 1 GeV2(blue long-dashed line) and = 10 GeV^ (red sohd line), and the AS 
integrand for = 1 GeV^(blue dot-dashed line) and = 10 GeV^ (red dashed line) integrands 
of Too for W ^90 GeV. 

is very sensitive, when varies between 1 GeV^ and 10 GeV^, to the overlapping of 
the dipole cross-section bandwidth (A > X^^^'{Q'^,W)) and the radial dipole distribution 
Pqq^\x); we are then probing with a high accuracy the quality of the shape of the dipole 
cross-section. 

At high the tails of the distributions plays a dominant role. The tail of the dis- 
tribution corresponds to the region where the integrand of the radial distribution can be 
approximated by an exponential fall. 



tail 



XKoiV^X) ~ > VXexp{-^X) , (6.22) 

where typically 

A-^i~— A_^«4.5. (6.23) 

In the case 

Q2 = 1 GeV^, the bandwidth of the dipole cross-section mostly overlaps 
with the peak of the distributions as A^'** (l,90) ~ ^peak^ while it only overlaps with the 
tail of the distribution when = ig QeV^, A^^* (10,90) ~ A*^'^ 

Figure 14 shows the normalized integrand of Too- This summarize our discussion on 
the respective role of the radial distribution and of the dipole cross-section. This integrand 
is peaked near the saturation radius r ~ 2Rq(x). Comparing the cases = 10 GeV^ and 

= 1 GeV^, we see that this peak is moving to the right as decreases, going through 
the bandwidth of the dipole cross-section. 

6.2 The radial distribution of the 7^ — )• px transition 

In figure 15 are shown respectively for = 1 GeV^ and = 10 GeV^, 

1. the full twist 3 (Total) rescaled radial distribution Pii(A, )U^(Q^)), where we distin- 
guish the two following contributions, 
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Figure 15: The Total Pii(A, A^MQ^)) (red solid lines) results and their WW (blue dot-dashed lines) 
and genuine (Gen) (orange dashed lines) contributions, as well as the AS (purple long-dashed line) 
result p['^^\\) and the dipole cross-sections a{x, A) (black dot-dot-dashed lines) at = 90 GeV'^, 
for Q2 ^ 1 GeV^ (thick lines) and = 10 GeV^ (thin lines). 



• the WW contribution P^f ^(A, /if,(Q2))^ 

• the genuine contribution Pf™(A, /i|,(Q^)), 

such as^ Pii(A, /x|.) = 

2. the asymptotic rescaled radial distribution P\^^\X), using the asymptotic distribu- 
tion amplitudes. 

Contrary to the 7|, — )• /Ol transition, we see that the dependence on n\ is quite strong if 
one compares Pii(A,/i|,) to p[i^\\). We cannot then restrict ourselves to only study the 
asymptotic case. 

It is interesting to estimate the average (A) obtained with the different distributions 
Pii, p[^^\ and -Pff , where p^^ and -Pf^" have been normalized separately. The 

explicit expression for the asymptotic distribution 

Pif\x) = ^V[f^ (A = ^ I dy[yyfl\y'+f)X'K,{V^yX), (6.24) 

leads to 

(A)Sf ) = j dAApii'(^^)(A) = ^ « 8.33. (6.25) 

®The tilde is to differentiate the contributions P^"^ to the distribution Pn from the distributions 
PjY^ or p^'p"' which are normalized separately. 



p^,^{\^^l) + p!^{\^l), 
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The average value of A estimated with the WW distribution is, 



JdXjdy {ywf^{v,^il)XK,{^X)) 
= g(48 + 13a|(;.^))^(A)if). (6.26) 

The effect of the term involving a}\{ii^p) in the r.h.s. of eq. (6.26) is under 4%, which 
indicates that P^^^X^ii^) ~ P{^{X,iiy) is a good approximation. The computation of 
the average values of A for the all the different contributions are given in table 4. 





Total 


WW 


genuine 


AS 


(A)ii (^^(1 GeV^)) 


6.3 


8.7 


3.2 


8.3 


(A),i (^.^(10 GeV^)) 


7.3 


8.5 


3.5 


8.3 



Table 4: Average values of (A) = (r Q) for the different contributions to the radial distribution for 
two values of fipiQ^)- 



The results in table 4 show that the 7^ — >■ pT transition is more sensitive to saturation 
effects than the 7^ — t- transition as (A)-,^^ is about twice larger than (A)qq. Indeed it 
means that more dipolcs arc produced in the bandwidth of the dipole cross-section by the 
radial distribution -Pii(A, p,^) than Poo(-^) A*!')- As a consequence the polarized cross-section 
or should be a more sensitive observable to probe features of the saturation regime than 

The transverse dipole scale associated to the WW contribution, using eq. (6.12), is 



Wii =Q\\fyy)i\ ~27' 



which is not so far from the values of the function pjp{Q ) = — ^ — - that is used here. 
The tail of the distribution Pj^^ {X, p'^p) can be defined as, 



A^ Ki ( V^A) ^ > A VX exp(- V^A) , (6.27) 



where 



^taiiw^ ^^^ ^10.9, (6.28) 

{vyy)ii 



in which the average value {y/yy)^^ is approximated by {y/yy) 11 ~ 0.37. 

The genuine contribution Pf^"(A, /i|.), which vanishes in the limit pj^ 00, depends 
strongly on the factorization scale p^{Q^), as one can see in figure 15. The distribution 
Pq"(A, /Up) can be split into two contributions, 

pr(A,Mi.) = p!r^''\x,f^i) + pfr^'''\x,f^i) , (6.29) 
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where 

Pfr^''\\t^l) = ^ / dyi,J^^'-'^\y,X/Q;Q,l^l), (6.30) 
P!r^'''\\f^h = ^ / dy£dy,i.;^^{-iy,,y,X/Q;Q,f^l). (6.31) 

pgen (gg) involves the exchange of two-parton in the hard part and the genuine so- 

lutions of the two-parton DAs, while Vfl^^ ^'^^^\x, luijp) involves the three-parton exchange 
hard part T'^^{yi,y2,L'iQ) defined in eq. (3.19). 

The transverse scale /j, present in Pf^"^^^^(A, /i|.) is of order '^^^^ (Mf) ~ 0.5 

and is not very sensitive to The choice of /Xp(Q^) ~ y is then a good choice for this 
contribution. 

The contribution Pf^^^'^^^\x, fxj^) have several transverse scales which are fii, fi2, f^qg, 
jiqg and defined by eqs. (3.20), each corresponding to a dipole configuration involv- 
ing two of the three partons available in the process. In order to estimate these trans- 
verse scales, we first evaluate the average fraction of momentum carried by the quark 
(yi)fr ^''"^ the antiquark (y2)fr ^''"'^ (/^l-) and the gluon (y,)fr ^^^'^^ (/x|), defined as 

I dXX [ dy2 I dyiyitpl^-^.'"^{yi,y2,r;Q,iJ,l) 

(y,)fr (/.I) = ^ \ ^° . (6.32) 

J dXX J dy2 dyiipJ^-~^.^'"^{yi,y2,r;Q,iJ,jr) 

Using the fact that (yi)fi'^*''^^^'* (m^) = {y2)Ti^'^^^^ il-'-l') due to the symmetry under the 
exchange of the quark and the antiquark, the obtained values are given in table 5. 



, > gen (qqg) , \ 


(2/2)11 [fJ'F) 




41.8% 


41.8% 


16.3% 



Table 5: Average values (yi)f™^'^''®' {iJ,jp) of the fraction of longitudinal momentum of the p meson 
carried by each of the parton i. 

An estimation of the transverse scales can be made, using 

^ = ^ ~ Vivi) (1 ' ivi)) « 0.49 , (6.33) 



(Nq) _ (/A 



Q Q 

ilMq) 



Q 



(yi) 








(yi) 


im) 







0.34 , (6.34) 
0.46 , (6.35) 



where we assume that the transverse scale values arc roughly approximated by using the 

average fractions of longitudinal momentum in eqs. (3.20). 

These values arc evaluated at ;u|^(l GeV^). Other values of fip have also been used, 

Q 

2 



leading to approximately the same results. We note that the function jJipiQ^) ~ % is 
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close to the {fii) and (fJ-qq) values, and of the same order of magnitude than (Hqy) but in 
principle, one should adapt the choice of the factorization scale to the relevant transverse 
scales at stake for each part of the process. 

The WW and the genuine contributions to the radial distribution Pii(A, are of the 
same order of magnitude for ~ 10 GeV^, and the genuine contribution becomes even 
more important for = 1 GeV^. At = 10 GeV^, as 

(A)ff (-3.5) < AS^*-(10GeV2)(~7.7) < (A)J^^ (~ 8.7) , 

most of the dipoles in the bandwidth of the dipole cross-section are provided by the WW 
contribution, which explains why the predictions for cjt arc dominated by the WW pre- 
dictions, as shown in figure 4(a), and consequently, why the results depend weakly on the 
factorization scale. 

However the fact that the genuine contribution is important even at large , indicates 
that the three-parton exchange between the and pT states is important in such relations 
as the normalization eq. (6.36) of the p meson wave function [7,81], and the electronic decay 
width eq. (6.37) [7,82], 

1 = E / ^^2/ y" |*5(y,r)|' , (6.36) 

h,h 

efpTUpie; ■ep) = ^jdyj (iV*^yy,r) ^^^(y,!) , (6.37) 
h,h 

where the exchange of only two-parton is assumed. Indeed, the r.h.s. of eq. (6.37), if one 
expands at large the p meson wave function around r = , is the WW approximated 
result, which therefore misses the genuine contributions arising from three-parton corre- 
lators, which can have a significant effect even for large values, see figure 15. These 

relations are usually used to constrain the parameters of the wave function models of the p 
meson, assuming that the meson is solely constituted of a quark and an antiquark, which 
then consist in neglecting the higher Fock state contributions like the genuine contribution. 
At = 1 GeV^, the peak of the genuine part of the radial distribution enters the dipole 
cross-section bandwidth (A^^* (l,90) < (A)!^"^) and gives an important contribution to the 
integrand of Tn. When increasing Q^, as shown in figure 16 where we display the product 
'Pii{i^:Q'^ , fJ''F)(^{x,r) for = 1 GeV^ and = 10 GcV^, we can note that the difference 
between the AS and the Total results is a consequence of the genuine contribution growing 
when decreases. It is the convolution with the dipole cross-section which washes-out 
the effect of these genuine twist-3 contributions. 

6.3 Comparison w^ith the radial distributions obtained from models of the p 
meson wave function. 

It is instructive to compare shapes of the radial distributions Vqo and Vu used in our 
analysis with those used in two other approaches which involve the overlap of the virtual 
photon wave functions and of the p meson wave functions: 
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Figure 16: The Total contributions at fJ.^{Q^) for = 1 GeV'^ (blue long-dashed line) and 
Q2 — iQ GeV^ (red soHd line), and the AS contributions for ~ 1 GeV^ (blue dot-dashed line) 
and = 10 GeV^ (red dashed line) to the normalized integrands of Tn, i.e. 'Pii{r, Q^, ^\)a{x, r), 
for W = QQ GeV. 



• the "Boosted Gaussian" (BG) model [46], 

• the"Gaus-LC" model [55]. 

We use here the convention and the parameter values of ref. [47], which for completeness 
are shown in table 6. The scalar parts of the wave functions are given by. 



Model 


Nt 


Rl GeV-2 


Nl 


Rl GeV-2 




Gaus-LC 


4.47 


21.9 


1.79 


10.4 


fp 


Boosted Gaussian 


0.911 


0.853 


12.9 




0.182 



Table 6: Parameter of the "Gaus-LG" and the "Boosted Gaussian" models taken from ref. [47], 
for Mp = 0.776 GeV, fp = 0.156 GeV, ruf = 0.14 GeV and with fj;^fp. 



„2 



^Gau..-LC(y^^) = Nt {yyfe ^ , (6.38) 

_ ^2 

^Gauss-LC(y^^)^^^y-g ^, (6.39) 

*l%(y.r) = N,,yy exp [-^ - ^ + . (6.40) 

The overlaps with the virtual photon wave function are, 

J^^'^y(y,r)M/J-(y,r:) oc m}Ko{fir)(t>T{y,r) - {y' + f)fiKi{f,r)drMy,r) , (6.41) 
h,h 

^^'^^-*(2/,r)^^^-(y,r:) « yyKo{fir) |mp</>L(y,r) + J> L(y, j , (6.42) 

h,h \ / 
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with 6 = for the Gaus-LC model and 6 = 1 for the BG model. The radial distributions 
thus read, 

rL,T{r) = j^rjdy^ (^^n) > (6-43) 

' h,h 

where the factors N'l,t normalize the distributions Vl,t{t~)- Comparing the distributions 




(c) 7y pT radial distributions and a at = (d) 7^ — > pr radial distributions and a at = 



1 GeV^ lOGeV^. 

Figure 17: The Gauss-LC (green, long dashed), BG (blue, dot-dot-dashed). Total (red, solid) 
and AS (red, dashed) radial distributions for the 7^ — > pL transition (top) and for the 7^ — >■ pT 
transition (bottom), vs r for = 1 GeV^ (left) and = 10 GeV^ (right), as well as the dipole 
cross-section a{x^r) rescaled by the factor 5(To for = 90 GeV (black, dot-dashed). 

at Q2 = 1 GeV2 (figures 17(a) and 17(c)) and = 10 GeV^ (figures 17(b) and 17(d)) 
we see that at large Q^^ in the bandwidth of the dipole cross-section, our distributions are 
converging with the distributions obtained from the Gaus-LC and BG models. 

In the 7^ — >• pt case in figures 17(c) and 17(d), we see that the BG and Gaus-LC models 
are closer to the distribution Vu than to the asymptotic distribution. When compared to 
other distributions, the asymptotic distribution V^^ is shifted to the right, thus selecting 
larger dipole sizes. 

In the 72 — )■ pL case, in figure 17(a) we see that the distributions from the Gaus-LC 
and BG models are not close to our predictions, indicating that the higher twist corrections 
are presumably more important at small than in the 7^ — )• pT transition. 

These qualitative remarks remains the same when using the AAMQS dipole models, 
and we expect that they are model independent. 
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7. Conclusions 

We performed phenomenological analysis of experimental data from HERA on p meson 
electroproduction within the approach based on the recently derived [44] impact parameter 
representation of 7* — t- p impact factor up to twist-3 accuracy. The important feature of 
this representation consists in the inclusion of contributions coming from both two- and 
three-partonic Fock states, maintaining a close connection with the dipole model picture. 
Consequently, it was possible to include in our framework the saturation effects. Our 
predictions show that we can get simultaneously good predictions for the polarized cross- 
sections (jt and gl. The abilility of the model to reproduce the data is the confirmation 
of the following points: 

• the factorization of the dipole cross-section in the helicity amplitudes of the elec- 
troproduction of the p meson works and, as an universal quantity, is the same for 
Too and Til, giving the good energy dependence and normalizations of the polarized 
cross-sections, 

• the collinear factorization procedure of the p meson is justified and works successfully 
beyond the leading twist. 

As expected the model has some limits due to the truncation of the twist expansion. 
Thanks to HERA data, we have identified the virtuality Q^™'" ~ 5 GeV^ were the higher 
twist corrections become important, which is a motivation to compute impact factors be- 
yond the twist 3 accuracy in order to probe the genuine saturation regime which starts at 
lGeV2. 

Other helicity amplitudes could be computed keeping the same approach, they would 
be useful in the t ^ tmin regime. The kinematics of the impact factor can be also extended 
to take into account the dependence of the impact factors, which would be a test for the 
dipole models which include the impact parameter dependence, providing a probe of the 
proton shape [7], in particular through local geometrical scaling [83,84]. 

Data also exist for leptoproduction. In this case quark-mass effects should be taken 
into account, in particular, because this allows the transversely polarized to couple 
through its chiral-odd twist-2 DA. Indeed, as it was pointed out in [42], the fact that 
the ratio Th/Toq is not the same (after trivial mass rescaling) for p and mesons points 
to the importance of this effect. This is also beyond the scope of our present study, but 
may open an interesting way for accessing chiral-odd DAs. 

The next-to-leading order effects - both on the evolution and on the impact factor - 
should be studied, since it is now known that both may have a important phenomenological 
effect [85-89]. 

On the experimental side, the future Electron-Ion Collider [90] and Large Hadron 
Electron Collider [91] with a high center-of-mass energy and high luminosities, as well as 
the International Linear Collider [92-94] will hopefully open the opportunity to study in 
more detail the hard diffractive production of mesons [85-87, 95-99] . 
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A. Appendices 

A.l Distribution amplitudes in the LCCF parametrization 

The seven chiral-even^ p-meson DAs up to twist 3 are defined through^ the following matrix 
elements of nonlocal light-cone operators [40] , for two-parton 

Jo 

{p{Pp) \4>{z)l^lpm\ 0> = trn,f,Rl^ [ dy e'yP-'ipA{y; p-l) , (A.2) 

Jo 

{p{Pp)\ i^izhf^i ^ ^-(0) |0) = mpfpp^e*T, [ dy e'^^'-^f (y; , (A.3) 

Jo 

{piPp)\ dl m |0) = i mpfpp^ Rl, [ dy e'^'-' ipl{y- ,,1) , (A.4) 

Jo 

and three-parton correlators 

(pipp) \^{zih^gAliz2)m\ 0) = mX(/^F)PM4a (A.5) 

X / d.y2 r dyie'y^P-'^+'^y^-y^^P-'^B{yi,y2;fil) 
Jo Jo 

(pip,) \^P{zih57Mli^2)m\ 0> =impfip{p,l)p^Rl„ (A.6) 

X [dy2 rdyie'y^P-'^+''^y'-y^^P-''D{y,,y2;p%), 
Jo Jo 

where we used the standard notation dp= ^{dp — dp) . 

DAs are linked by linear differential relations derived from equations of motion and 
n— independency condition [39,40]. The solutions for ^pp{y) = {ip^,ipA-,^'i ,v'a} 
sum of the solutions in the so-called WW approximation and of genuine solutions, 

V.p(y) = v^ir^(y) + ^r(y)- (A.7) 

'^The chiral-odd twist-2 DA for the transversely polarized p meson does not contribute to the process 
considered at the accuracy discussed here. This is also true in the approach based on coUinear factorization 
of generalized parton distributions [100,101]. 

®In the approximation where the mass of the quarks is neglected with respect to the mass of the p meson. 
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The WW approximation consists in neglecting the contribution from three-parton opera- 
tors, thus taking B{yi,y2] = D{yi,y2; fijp) = 0. Then, (fi^^ (y) become functions of 
ipi{y) only, and their explicit expressions are given by 



^2. {y^fJ'F) = 



u 



Jy 

I 

" Jy 



u 



u 



y l\u^^^^^ + y f\u^^^^^ 

U Jy U 



Jo " Jy 



u 



(A.8) 
(A.9) 
(A.IO) 
(A.ll) 



Genuine solutions only depend on {B{yi,y2; jJ^p), D{yi,y2] /^|')} or equivalently on the com- 
binations {5(yi,y2;MF))-^(yi)y2;MF)} defined by eq. (3.3), namely 



¥'f"(y;AiF) = ^ 



du h / du- 

y Jy 



i 

Jy U Jy 



U 

■ „2 



du 



A(u; fij,) 



u 



(A.12) 
(A.13) 



where A(u] has the compact form 



f'U 

A(u; lip) = / dy2 
Jo 



1 



9„ 



.y2 - u 

and it obeys the conditions 



/■I 1 

M{y2,u;iil)+ dy2 M{u,y2; fij.) (A.14) 

Ju y2-u 



I du A{u; lip) = and [ duuA{u;np) = 0, (A.15) 
Jo Jo 



coming, respectively, from the constraints 



/ '/'fly; = and / {y - y) if^^'^iy; nj,) dy = . (A. 16) 

Jo Jo 

Equations (A.12) and (A.13) determine the expressions of ipf^ {y, fJ'^) and ip^^^{y; li^) as 



fi^^''{y;i4') = / duipf''{u;ii]r) 
Jo 



1 f\y^ dy /^^'' ^ ^(^1.^2; III) 
^ Jy ^ y2-yi 



T gen 

"Pa 



ry 

(v^I^f) = / dU'fA"'{u;iijp) 
Jo 



(A.17) 



(A.18) 



-ll^dy^l dy2 



\ S{yi,y2;iJ'l) - M{yi,y2;iPp) 



y2 - yi 
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The correspondence between our set of DAs and the one defined in ref. [50] is achieved 
through the fohowing dictionary derived in ref. [40]. It reads, for the two-parton vector 
DAs, 



and for the axial DA, 



1 dg^^\y) 



VAiy) 



dy 



For the three-parton DAs, the identification is 



B{y,,y.) = -'^^^^^^ and D^y,, y.) = -^^^^ ' ' 



(A.19) 



(A.20) 



(A.21) 



y2 -yi y2- yi 

ExpUcit forms for tpi, B, and D are obtained with the help of the results of ref. [50] obtained 
within the QCD sum rules approach. The first terms of the expansion in the momentum 
fractions of the three independent DAs thus have the form 

3, 



^i{y,i^F) = ^yy 



l + aliiil)-i5{y-yf-l) 



Biyi,y2;fJ'F) = -5040yiy2(yi - y2)iy2 - yi] 



Diyi,y2;fJ'F) = -360?/iy2(2/2 - yi) 



CO 



1 + 



{1,0} 



(7(y2 - yi) - 3) 



(A.22) 
(A.23) 
(A.24) 



The dependences on the renormalization scale /x^ of the coupling constants a^, u)^^ q|, C^, 
and ^i'^ given in ref. [50]. In appendix A.2 we present both the evolution equations 
and the values of these constants at = 1 GeV^ used in our analysis, as well as the 
dependence on fxp oi the DAs. 

A.2 Evolutions of DAs and coupling constants with the renormalization scale 

The parameters entering the DAs at = 1 GeV^ are updated^ in ref. [51] and their 
evolution equations are given in ref. [50] , we recall in table 7 their values for the p meson. 
For al, the evolution equation is 

a|(/x2) = al(/x2)L(/.2p/^o 



with 



where bo = {IIN^ - 2Nf)/3 , 7„ = 4Cf (^tp{n + 2) + 7^ - ^ - - 



1 



2(n+l)(n+2) 



— 7b + Ylk=i ^/k ■ For the f^p coupling constant, the evolution is given by 



(A.25) 

(A.26) 
and ip{n) = 

(A.27) 



®We use the notations of ref. [50] for the parameters, they are related to the updated parameters of 
ref. [51] by the following relations, C3* = Cl, Cl' = '*^3/14 and Cl'^fi.o} = '^3- 
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0.52 


'^{1,0} 


-3.0 


1 y 


28/3 


a| 


0.15 


S3 


0.030 


S3 


0.011 



Table 7: Coupling constants and Gegenbauer coefficients entering the p meson DAs, at the scale 
/Uo = 1 GeV updated in ref. [51]. Note that in ref. [50] the normalization are such that j^q^ — 

rV,A 

"^P 3p [here] ' 



with To = 



evolution equation [50]. Defining 



3Cg {Cg = Nc). The couplings f^p and 0;^^ i}(a*^)/3^(a*^) enter a matrix 



it reads 



with given by 



J- q 



''-K,l]/3^p(/^')-'-fo,l}('"')/3l,(M^)' 
^ |Cf + \Cg \Cf - \Cg ^ 



(A.28) 



(A.29) 



(A.30) 



Hence we get the dependence of /g^ and a;|J) ^ by diagonalizing the system. The dependence 
of the DAs on the renormalization scale is shown in ref. [42]. The DAs exhibit a non- 
negligible effect of QCD evolution, in particular, for the genuine twist-3 contributions. We 
recall that we chose the coUinear factorization scale of production of the p meson /Xi? to be 
equal to the renormalization scale of the process thus the dependence of the coupling 
constant in np is given by eqs. (A.25, A.27, A.29). 

A. 3 Dipole-proton scattering amplitude in the GS-Model 

In this appendix we calculate the Fourier transform of the proton impact factor. We denote 
by r = |r| the transverse size of the dipole. We should thus compute 



M2 + /fc2 



—ik-r\ 



(A.31) 
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Using the notation k = \k\ , the angular integration leads to 

a{r) = AnAj dk^^ - j^^) (1 - Mkr)) . (A.32) 



Relying on the identity 

1 /I k 



P M'^ + k'^J W^yk M^ + k'^J' 

we rewrite a{r) as 



(A.33) 



<^ir) = -j^I (A.34) 

with 

= /i - /2 - Is + /4 . (A.35) 

The integral I4 is UV and IR finite and reads [102] 

l4 = Ko{Mr). (A.36) 

The integrals Ii and I2 are both IR divergent and are regularized through dimensional 
regularization, while the UV divergencies of Ii and I3 are regularized through a cut-off A . 
We thus write 

A = /^^-^ + lnA. (A.37) 

Using the relation [102] 

r Ma x) dx = 2^a^^-^ + + i^) (^.38) 

Jo r (2 + 2^^-2^) 

we obtain 

^2 = y^ ^Jo(A;r) = 2-i+V-^;^^^~--lnr-7 + ln2. (A.39) 



Finally, 



^3 = [ Jr^'^'^ ~ In A - InM . (A.40) 



Thus, combining eqs. (A.37, A.39, A.40, A.36) we get 



M r 

/= (7 + ln— + iCo(Mr) ) , (A.41) 



and thus 

iirA/ , Mr 



<^{r) = ^[l + lT^^ + Ko{Mr)] . (A.42) 
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A.4 Results using the GBW and AAMQSb models 

We present some of the predictions obtained by using the GBW or the AAMQSb model 
for the dipole cross-section. As expected the results are not so far from the ones obtained 
with the AAMQSa model. In figure 18 and 19 are respectively shown for the AAMQSb 
and the GBW model, the polarized cross-sections or and gt- The spin density matrix 
element rpo predictions using these dipole models are shown in figure 20, for completeness 
we show also the prediction obtained with the GS model. 
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(7r(nb) 




(a) (TT vs HI data [4]. 



c^L(nb) 




(b) (TL VS HI data [4]. 



Figure 18: AS (purple dashed lines), WW (blue long dashed lines) and Total (red solid lines) 
contributions to ctt and ctl vs Q^, for W = 75 GeV, using the AAMQSb-model, compared to the 
data of HI [4]. 
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(7r(nb) 




(a) (JT vs HI data [4]. 



c^L(nb) 




(b) (TL VS HI data [4]. 



Figure 19: AS (purple dashed lines), WW (blue long dashed lines) and Total (red solid lines) 
contributions to ctt and ci, vs Q^, for W = 75 GeV, using the GBW-model, compared to the data 
of HI [4]. 
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(c) Too contributions using GBW-model vs HI data. 
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(d) r-Qo contributions using GBW-model vs ZEUS 
data. 
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Figure 20: Predictions for Tqq vs W and compared respectively with HI [4] and ZEUS [3] data, 
the AS (purple dashed lines), WW (blue long dashed lines). Total (red solid lines) contributions 
are shown separately using the AAMQSb-model, the GBW-model or the GS-model. 



References 

[1] ZEUS Collaboration, J. Breitweg et. al., Exclusive electroproduction of and J/ip mesons 
at HERA, Eur. Phys. J. C6 (1999) 603-627, [hep-ex/9808020]. 

[2] ZEUS Collaboration, J. Breitweg et. al, Measurement of the spin-density matrix elements 
in exclusive electroproduction of mesons at HERA, Eur. Phys. J. C12 (2000) 393-410, 
[hep-ex/9908026]. 

[3] ZEUS Collaboration, S. Chekanov et. al., Exclusive p° production in deep inelastic 
scattering at HERA, PMC Phys. Al (2007) 6, [arXiv : 0708 . 1478]. 



41 



[4] HI Collaboration, F. D. Aaron et. at, Diffractive Electroproduction of p and (j) Mesons at 
HERA, JHEP 05 (2010) 032, [arXiv : 0910 . 5831]. 

[5] A. Donnachie and P. Landshoff, Exclusive p Production in Deep Inelastic Scattering, 
Phys.Lett. B185 (1987) 403. 

[6] J. Nemchik, N. N. Nikolaev, E. Predazzi, and B. Zakharov, Color dipole systematics of 
diffractive photoproduction and electroproduction of vector mesons, Phys.Lett. B374 (1996) 
199-204, [liep-pli/9604419]. 

[7] S. Munier, A. M. Stasto, and A. H. Mueller, Impact parameter dependent S-matrix for 
dipole proton scattering from, diffractive meson electroproduction, Nucl. Phys. B603 (2001) 
427 445, [hep-ph/0102291]. 

[8] H. Cheng and T. T. Wu, Photon-photon scattering close to the forward direction, Phys. Rev. 
Dl (1970) 3414-3415. 

[9] G. Frolov and L. Lipatov Sov. J. Nucl. Phys. 13 (1971) 333. 

[10] V. Gribov, L. Lipatov, and G. Frolov, The leading singularity in the j plane in quantum 
electrodynamics, Sov.J.Nucl.Phys. 12 (1971) 543. 

[11] S. Catani, M. Ciafaloni, and F. Hautmann, Gluon contributions to small x heavy flavor 
production, Phys. Lett. B242 (1990) 97. 

[12] S. Catani, M. Ciafaloni, and F. Hautmann, High-energy factorization and small x heavy 
flavor production, Nucl. Phys. B366 (1991) 135-188. 

[13] J. C. Collins and R. K. Ellis, Heavy quark production in very high-energy hadron collisions, 
Nucl. Phys. B360 (1991) 3-30. 

[14] E. M. Lcivin. M. G. Ryskin, Y. M. Shabclski, and A. G. Shuvaev, Heavy quark production in 
semihard nucleon interactions, Sov. J. Nucl. Phys. 53 (1991) 657. 

[15] V. S. Fadin, E. A. Kuraev, and L. N. Lipatov, On the Pomeranchuk Singularity in 
Asymptotically Free Theories, Phys. Lett. B60 (1975) 50-52. 

[16] E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, Multi - Reggeon Processes in the Yang-Mills 
Theory, Sov. Phys. JETP 44 (1976) 443-450. 

[17] E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, The Pomeranchuk Singularity in Nonabelian 
Gauge Theories, Sov. Phys. JETP 45 (1977) 199-204. 

[18] I. I. Balitsky and L. N. Lipatov, The Pomeranchuk Singularity in Quantum 
Chromodynamics, Sov. J. Nucl. Phys. 28 (1978) 822-829. 

[19] V. S. Fadin, R. Fiore, and M. L Kotsky, Gluon Regge trajectory in the two-loop 
approximation, Phys. Lett. B387 (1996) 593-602, [hep-ph/9605357]. 

[20] G. Camici and M. Ciafaloni, Irreducible part of the next-to-leading BFKL kernel, Phys. Lett. 
B412 (1997) 396-406, [hep-ph/9707390]. 

[21] M. Ciafaloni and G. Camici, Energy scale(s) and next-to-leading BFKL equation, Phys. Lett. 
B430 (1998) 349-354, [hep-ph/9803389]. 

[22] V. S. Fadin and L. N. Lipatov, BFKL pom.eron in the next-to-leading approximation, Phys. 
Lett. B429 (1998) 127-134, [hep-ph/9802290]. 



42 



[23] A. H. Mueller, Small x Behavior and Parton Saturation: A QCD Model, Nucl. Phys. B335 
(1990) 115. 

[24] N. N. Nikolaev and B. G. Zakharov, Colour transparency and scaling properties of nuclear 
shadowing in deep inelastic scattering, Z. Phys. C49 (1991) 607-618. 

[25] S. J. Brodsky, L. Frankfurt, J. F. Gunion, A. H. Mueller, and M. Strikman, Diffractive 
leptoproduction of vector mesons in QCD, Phys. Rev. D50 (1994) 3134-3144, 
[hep-ph/9402283]. 

[26] L. Frankfurt, W. Koepf, and M. Strikman, Hard diffractive electroproduction of vector 
m.esons in QCD, Phys. Rev. D54 (1996) 3194 3215, [hep-ph/9509311]. 

[27] J. C. Collins, L. Frankfurt, and M. Strikman, Factorization for hard exclusive 

electroproduction of mesons in QCD, Phys. Rev. D56 (1997) 2982-3006, [hep-ph/96 11433]. 

[28] A. V. Radyushkin, Nonforward parton distributions, Phys. Rev. D56 (1997) 5524-5557, 
[hep-ph/9704207]. 

[29] G. R. Farrar and D. R. Jackson, The Pion Form-Factor, Phys. Rev. Lett. 43 (1979) 246. 

[30] G. P. Lepage and S. J. Brodsky, Exclusive Processes in Quantum Chromodynamics: 

Evolution Equations for Hadronic Wave Functions and the Form-Factors of Mesons, Phys. 
Lett. B87 (1979) 359-365. 

[31] A. V. Efremov and A. V. Radyushkin, Factorization and Asymptotical Behavior of Pion 
Form-Factor in QCD, Phys. Lett. B94 (1980) 245-250. 

[32] D. Yu. Ivanov, L. Szymanowski, and G. Krasnikov, Vector meson electroproduction at 
next-to-leading order, JETP Lett. 80 (2004) 226-230, [hep-ph/0407207]. 

[33] H.-N. Li and G. Sterman, The Perturbative pion form-factor with Sudakov suppression, 
Nucl. Phys. B381 (1992) 129-140. 

[34] M. Vanderhaeghen, P. A. M. Guichon, and M. Guidal, Deeply virtual electroproduction of 
photons and mesons on the nucleon: Leading order amplitudes and power corrections, Phys. 
Rev. D60 (1999) 094017, [hep-ph/9905372]. 

[35] S. V. Goloskokov and P. Kroll, Vector meson electroproduction at small Bjorken-x and 
generalized parton distributions, Eur. Phys. J. C42 (2005) 281-301, [hep-ph/0501242]. 

[36] S. V. Goloskokov and P. Kroll, The longitudinal cross section of vector meson 
electroproduction, Eur. Phys. J. C50 (2007) 829-842, [hep-ph/06 11290]. 

[37] S. V. Goloskokov and P. Kroll, The role of the quark and gluon GPDs in hard vector-meson 
electroproduction, Eur. Phys. J. C53 (2008) 367-384, [/08.3569]. 

[38] D. Yu. Ivanov and R. Kirschner, Polarization in diffractive electroproduction of light vector 
mesons, Phys. Rev. D58 (1998) 114026, [hep-ph/9807324]. 

[39] I. V. Anikin, D. Yu. Ivanov, B. Pire, L. Szymanowski, and S. Wallon, On the description of 
exclusive processes beyond the leading twist approximation, Phys. Lett. B682 (2010) 
413-418, [arXiv: 0903. 4797]. 

[40] I. V. Anikin, D. Yu. Ivanov, B. Pire, L. Szymanowski, and S. Wallon, QCD factorization of 
exclusive processes beyond leading twist: 7^ — > pr impact factor with twist three accuracy, 
Nucl. Phys. B828 (2010) 1-68, [arXiv: 0909. 4090]. 



43 



[41] I. F. Ginzburg, S. L. Panfil, and V. G. Scrbo, Possibility of the experimental investigation of 
the QCD Pomeron in semihard processes at the gamma gamma collisions, Nucl. Phys. B284 
(1987) 685-705. 

[42] I. V. Anikin, A. Besse, D. Yu. Ivanov, B. Pire, L. Szymanowski, and S. Wallon, A 

phenomenological study of helicity amplitudes of high energy exclusive leptoproduction of the 
p meson, Phys. Rev. D84 (2011) 054004, [arXiv: 1105. 1761]. 

[43] J. F. Gunion and D. E. Soper, Quark Counting and Hadron Size Effects for Total Cross- 
Sections, Phys. Rev. D15 (1977) 2617 2621. 

[44] A. Besse, L. Szymanowski, and S. Wallon, The dipole representation of vector meson 

electroproduction beyond leading twist, Nucl. Phys. B867 (2013) 19-60, [arXiv: 1204.2281]. 

[45] J. Nemchik, N. N. Nikolaev, E. Predazzi, and B. G. Zakharov, Color dipole phenomenology 
of diffractive electroproduction of light vector mesons at HERA, Z. Phys. C75 (1997) 71-87, 
[hep-ph/9605231]. 

[46] J. R. Forshaw, R. Sandapen, and G. Shaw, Colour dipoles and p, $ electroproduction, Phys. 
Rev. D69 (2004) 094013, [hep-ph/0312172]. 

[47] H. Kowalski, L. Motyka, and G. Watt, Exclusive diffractive processes at HERA within the 
dipole picture, Phys. Rev. D74 (2006) 074016, [hep-ph/0606272]. 

[48] J. R. Forshaw and R. Sandapen. Extracting the p meson wavefunction from HERA data, 
JHEP 11 (2010) 037, [arXiv: 1007. 1990]. 

[49] J. Forshaw and R. Sandapen, Extracting the Distribution Amplitudes of the p meson from 
the Color Glass Condensate, JHEP 1110 (2011) 093, [arXiv: 1104.4753]. 

[50] P. Ball, V. M. Braun, Y. Koike, and K. Tanaka, Higher twist distribution amplitudes of 
vector mesons in QCD: Formalism and twist three distributions, Nucl. Phys. B529 (1998) 
323-382, [hep-ph/9802299]. 

[51] P. Ball, V. Braun, and A. Lcnz, Twist-4 distribution amplitudes of the K* and <p mesons in 
QCD, JHEP 0708 (2007) 090, [arXiv:0707. 1201]. 

[52] K. J. Golcc-Biernat and M. Wusthoff, Saturation effects in deep inelastic scattering at low 
Q"^ and its implications on diffraction, Phys. Rev. D59 (1999) 014017, [hep-pli/9807513]. 

[53] J. L. Albacete and Y. V. Kovchegov, Solving high energy evolution equation including 
running coupling corrections, Phys. Rev. D75 (2007) 125021, [arXiv: 0704. 0612]. 

[54] J. L. Albacete, Particle multiplicities in Lead-Lead collisions at the LHC from non-linear 
evolution with running coupling, Phys. Rev. Lett. 99 (2007) 262301, [arXiv: 0707. 2545]. 

[55] H. Kowalski and D. Teaney, An Impact parameter dipole saturation model, Phys. Rev. D68 
(2003) 114005, [hep-ph/0304189]. 

[56] J. R. Forshaw and D. A. Ross, Quantum chromodynamics and the pomeron, Cambridge 
Led. Notes Phys. 9 (1997) 1-248. 

[57] G. Frolov, V. Gribov, and L. Lipatov, On Regge poles in quantum electrodynamics, 
Phys. Lett. B31 (1970) 34. 

[58] V. S. Fadin and A. D. Martin, Infrared safety of impact factors for colorless particle 
interactions, Phys. Rev. D60 (1999) 114008, [hep-ph/9904505]. 



44 



[59] The HERMES Collaboration, A. Airapctian et. at, Ratios of Helicity Amplitudes for 
Exclusive fp Electroproduction, Eur. Phys. J. C71 (2011) 1609, [arXiv : 1012 . 3676]. 

[60] A. G. Shuvaev, K. J. Golec-Biernat, A. D. Martin, and M. G. Ryskin, Off-diagonal 

distributions fixed by diagonal partons at small x and ^, Phys. Rev. D60 (1999) 014015, 
[hep-ph/9902410]. 

[61] A. D. Martin, M. G. Ryskin, and T. Teubncr, dependence of diffractive vector meson 
electroproduction, Phys. Rev. D62 (2000) 014022, [hep-ph/9912551]. 

[62] K. J. Golec-Biernat and M. Wusthoff, Saturation in diffractive deep inelastic scattering, 
Phys. Rev. D60 (1999) 114023, [hep-ph/9903358]. 

[63] E. Gotsman, E. Levin, M. Lublinsky, and U. Maor, Towards a new global QCD analysis: 
Low X DIS data from nonlinear evolution, Eur. Phys. J. C27 (2003) 411-425, 
[hep-ph/0209074]. 

[64] J. Bartcls, K. J. Golcc-Bicrnat, and H. Kowalski, A modification of the saturation model: 
DGLAP evolution, Phys. Rev. D66 (2002) 014001, [hep-ph/0203258]. 

[65] E. lancu, K. Itakura, and S. Munier, Saturation and BFKL dynamics in the HERA data at 
small X, Phys. Lett. B590 (2004) 199-208, [hep-ph/0310338]. 

[66] J. L. Albacete, N. Armesto, J. G. Milhano, C. A. Salgado, and U. A. Wiedemann, Nuclear 
size and rapidity dependence of the saturation scale from QCD evolution and experimental 
data, Eur.Phys.J. C43 (2005) 353-360, [hep-ph/0502167]. 

[67] V. Goncalves, M. Kugeratski, M. Machado, and F. Navarra, Saturation physics at HERA 
and RHIC: An Unified description, Phys.Lett. B643 (2006) 273- 278, [hep-ph/0608063]. 

[68] J. Jalilian-Marian, A. Kovner, A. Leonidov, and H. Weigert, The Wilson renormalization 
group for low x physics: Towards the high density regime, Phys. Rev. D59 (1999) 014014, 
[hep-ph/9706377]. 

[69] J. Jalilian-Marian, A. Kovner, and H. Weigert, The Wilson renormalization group for low x 
physics: Gluon evolution at finite parton density, Phys. Rev. D59 (1999) 014015, 
[hep-ph/9709432]. 

[70] A. Kovner, J. G. Milhano, and H. Weigert, Relating different approaches to nonlinear QCD 
evolution at finite gluon density, Phys. Rev. D62 (2000) 114005, [hep-ph/0004014]. 

[71] H. Weigert, Unitarity at small Bjorken x, Nucl. Phys. A703 (2002) 823-860, 
[hep-ph/0004044]. 

[72] E. lancu, A. Leonidov, and L. D. McLcrran, Nonlinear gluon evolution in the color glass 
condensate. I, Nucl. Phys. A692 (2001) 583-645, [hep-ph/0011241]. 

[73] E. Ferreiro, E. lancu, A. Leonidov, and L. McLcrran, Nonlinear gluon evolution in the color 
glass condensate. H, Nucl. Phys. A703 (2002) 489- 538, [liep-pli/0109115]. 

[74] I. Balitsky, Operator expansion for high-energy scattering, Nucl. Phys. B463 (1996) 99-160, 
[liep-pli/9509348]. 

[75] Y. V. Kovchegov, Small-x F2 structure function of a nucleus including multiple pomeron 
exchanges, Phys. Rev. D60 (1999) 034008, [hep-ph/990128l]. 

[76] J. L. Albacete, N. Armesto, J. G. Milhano, and C. A. Salgado, Non-linear QCD meets data: 
A Global analysis of lepton-proton scattering with running coupling BK evolution, Phys. Rev. 
D80 (2009) 034031, [arXiv: 0902. 1112]. 



45 



[77] J. L. Albacctc, N. Armcsto, J. G. Milhano, P. Quiroga-Arias, and C. A. Salgaclo, AAMQS: 
A non-linear QCD analysis of new HERA data at small-x including heavy quarks^ 
Eur.Phys.J. C71 (2011) 1705, [arXiv: 1012.4408]. 

[78] L. D. McLerran and R. Venugopalan, Boost covariant gluon distributions in large nuclei, 
Phys.Lett. B424 (1998) 15-24, [nucl-th/9705055]. 

[79] A. H. Mueller, Soft gluons in the infinite momentum wave function and the BFKL pomeron, 
Nucl. Phys. B415 (1994) 373-385. 

[80] H. Navelet, R. B. Pcschanski, C. Royon, and S. Wallon, Proton structure functions in the 
dipole picture of BFKL dynamics, Phys. Lett. B385 (1996) 357-364, [hep-ph/9605389]. 

[81] G. P. Lepage and S. J. Brodsky, Exclusive Processes in Perturbative Quantum 
Chromodynamics, Phys. Rev. D22 (1980) 2157. 

[82] H. G. Dosch, T. Gousset, G. Kulzinger, and H. J. Pirner, Vector meson leptoproduction and 
nonperturbative gluon fluctuations in QCD, Phys. Rev. D55 (1997) 2602-2615, 
[liep-pli/9608203]. 

[83] E. Ferreiro, E. lancu, K. Itakura, and L. McLerran, Froissart bound from gluon saturation, 
Nucl.Phys. A710 (2002) 373-414, [hep-ph/020624l]. 

[84] S. Munier and S. Wallon, Geometric scaling in exclusive processes, Eur. Phys. J. C30 
(2003) 359-365, [hep-ph/0303211]. 

[85] D. Yu. Ivanov and A. Papa, Electroproduction of two light vector mesons in the next- 
to-leading approximation, Nucl. Phys. B732 (2006) 183-199, [hep-ph/0508162]. 

[86] D. Yu. Ivanov and A. Papa, Electroproduction of two light vector mesons in next-to-leading 
BFKL: Study of systematic effects, Eur. Phys. J. C49 (2007) 947-955, [liep-ph/06 10042]. 

[87] F. Caporale, A. Papa, and A. S. Vera, Collinear improvement of the BFKL kernel in the 
electroproduction of two light vector mesons, Eur. Phys. J. CSS (2008) 525-532, 
[arXiv: 0807. 0525]. 

[88] D. Colferai, F. Schwennsen, L. Szymanowski, and S. Wallon, Mueller Navelet jets at LHC - 
complete NLL BFKL calculation, JHEP 12 (2010) 026, [arXiv: 1002. 1365]. 

[89] B. Ducloue, L. Szymanowski, and S. Wallon, Mueller- Navelet jets at LHC: the first complete 
NLL BFKL study, PoS" QNP2012 (2012) 165, [arXiv: 1208.6111]. 

[90] D. Boer, M. Diehl, R. Milner, R. Venugopalan, W. Vogelsang, et. ai, Gluons and the quark 
sea at high energies: Distributions, polarization, tomography, arXiv: 1108. 1713. 

[91] LHeC Study Group Collaboration, J. Abelleira Fernandez et. ai, A Large Hadron 

Electron Collider at CERN: Report on the Physics and Design Concepts for Machine and 
Detector, J.Phys. G39 (2012) 075001, [arXiv: 1206.2913]. 

[92] ILC Collaboration, J. Brau, (Ed. ) et. ai, LLC Reference Design Report Volume 1 - 
Executive Summary, arXiv: 0712. 1950. 

[93] ILC Collaboration, G. Aarons et. ai. International Linear Collider Reference Design 
Report Volume 2: PHYSICS AT THE ILC, arXiv : 0709 . 1893. 

[94] ILC Collaboration, T. Behnke, (Ed. ) et. al., ILC Reference Design Report Volume 4 - 
Detectors, arXiv : 0712 . 2356. 



46 



[95] B. Pire, L. Szymanowski, and S. Wallon, Double diffractive p-production in 7*7* collisions, 
Eur. Phys. J. C44 (2005) 545 558, [liep-pli/0507038]. 

[96] R. Enberg, B. Pire, L. Szymanowski, and S. Wallon, BFKL resummation effects in 
7*7* pp, Eur. Phys. J. C45 (2006) 759-769, [hep-ph/0508134]. 

[97] B. Pire, M. Segond, L. Szymanowski, and S. Wallon, Q CD factorizations in 7*7* PlPl^ 
Phys. Lett. B639 (2006) 642-651, [hep-ph/0605320]. 

[98] M. Segond, L. Szymanowski, and S. Wallon, Diffractive production of two p^ mesons in 
e+e- collisions, Eur. Phys. J. C52 (2007) 93-112, [hep-ph/0703166]. 

[99] M. Segond, L. Szymanowski, and S. Wallon, A test of the BFKL resummation at ILC. Acta 
Phys. Polon. B39 (2008) 2577-2582, in Proceedings of the School on QCD, low-x physics, 
saturation and diffraction, Copanello Calabria, Italy, July 1-14 2007, [arXiv: 0802. 4128]. 

[100] M. Diehl, T. Gousset, and B. Pire, Exclusive electroproduction of vector mesons and 
transversity distributions, Phys. Rev. D59 (1999) 034023, [hep-ph/9808479]. 

[101] J. C. Collins and M. Diehl, Transversity distribution does not contribute to hard exclusive 
electroproduction of mesons, Phys. Rev. D61 (2000) 114015, [hep-ph/9907498] . 

[102] I. Gradshtcyn and I. Ryzhik, Table of Integrals, Series and Products, Academic Press, New 
York (1980). 



47 



